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and  electron  response  to  channel  D1  (middle  and  right  panels). 

2.  Shown  here  (left  to  right)  in  a  format  identical  to  Figure  1  is  the  proton  response  for 
channels  D02  and  T02,  and  the  electron  response  for  T06. 
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1.  INTRODUCTION 


A  previously  published  technical  report  [Brautigam,  et  al.,  2006;  hereafter  referred  to  as  TR#i  ] 
was  the  first  in  a  series  by  AFRL  to  characterize  the  response  of  the  Compact  Environmental 
Anomaly  Sensor  (CEASE).  That  report  (TR#1 )  includes  (1 )  a  definition  of  the  response  function 
and  its  relation  to  the  geometric  factor;  (2)  a  description  of  the  Monte  Carlo  simulations  that  were 
performed  to  determine  the  proton  and  electron  response  functions  for  the  telescope  and  dosimeters; 
and  (3)  the  resulting  response  functions  and  a  comparison  with  experimental  calibration  results. 

This  report  details  how  these  response  functions  are  used  to  estimate  a  set  of  geometric  factors  for  a 
set  of  integral  flux  channels  derived  for  both  the  telescope  and  dosimeters. 

As  described  in  TR#1,  CEASE  includes  a  telescope  and  two  dosimeters  (thin  shield,  DDl;  thick 
shield,  DD2).  The  telescope  includes  a  front  (DFT)  and  back  (DBT)  solid  state  detector.  A  particle 
passing  through  the  telescope  will  deposit  some  portion  of  its  incident  energy  (either  0  or  some  non¬ 
zero  fraction)  in  each  of  the  2  detectors.  The  energy  deposition  scale  for  each  of  DFT  and  DBT  is 
divided  into  9  intervals  so  that  the  pair  of  detectors  leads  to  81  (9x9)  different  logic  blocks 
designated  by  LB(1,J),  where  1  and  J  label  the  energy  deposition  interval  for  DBT  and  DFT, 
respectively.  Various  groups  of  adjacent  logic  blocks  are  summed  on  board  the  satellite  to  reduce 
telemetry  demands.  Each  dosimeter  (DDl  and  DD2)  has  three  Linear  Energy  Transfer  (LET) 
channels  defined  by  the  energy  deposition  range  in  the  dosimeter's  single  solid  state  detector,  which 
for  CEASE  sn002,  are  the  following:  LoLET  (0.05  -  0.8  MeV),  HiLETA  (0.8  -  3.0  MeV)  and 
HiLETB  (3.0-1 0.  MeV),  for  both  flux  count  rates  and  dose  count  rates.  For  CEASE  sn007,  the 
energy  deposition  boundary  between  LoLET  and  HiLETA  channels  is  1.0  MeV  instead  of  0.8  MeV. 
The  telescope  logic  block  channels  (individual  or  summed)  and  the  dosimeter  channels  that  are 
included  in  the  telemetry  stream  are  listed  in  Appendix  A. 

2.  CEASE  GEOMETRIC  FACTORS 

An  isotropic  differential  flux  spectrum  j(E)  is  related  to  the  i”’  channel  count  rate  Ci  according  to 
the  following  integral,  where  Rj  (E)  is  the  energy  dependent  response  function  for  the  i*'’  channel; 

C,=\RXE)j{E)dE  (1) 

0 

A  consistent  set  of  units  for  the  above  variables  are:  Cj  [counts/s],  Ri(E)  [cm^  sr],  E  [MeV],  and  j(E) 
[#/(cm^  s  sr  MeV)].  In  the  case  where  Rj(E)  can  be  approximated  by  a  step  function  with  an  energy 
threshold  Et  and  constant  non-zero  value  R|(E)  =  Gj  for  E  >  Et,  (1 )  can  be  simplified  to: 

rn 

C,=G,jj(E)dE 

With  the  definition  of  integral  flux  =  |  j{E)dE ,  Eq.  (2)  reduces  to: 


(2) 


In  the  case  where  Ri(E)  can  be  approximated  by  a  constant  non-zero  value  (g,)  within  a  relatively 
narrow  energy  interval  [£■/,  £2^  throughout  which  j(E)  - 7,  =  constant,  then  Eq.  (1)  reduces  to 
C;  =  g,  and  the  differential  flux  for  the  /*  channel  [with  E=(Ei+E2)/2  ;  AE,  =  E2-  Ei]  is 
given  by: 


J>iE)  = 


C, 


(4) 


where j(E)  is  a  differential  flux  in  units  of  (cm“  s  sr  MeV)  ’. 

In  principle,  given  a  set  of  channel  count  rates  and  response  functions,  one  can  invert  the  set  of 
corresponding  integral  equations  Eq.  (1)  to  solve  for  the  unknown  differential  flux.  However,  in 
practice  this  is  problematic  because  one  typically  is  faced  with  an  “ill-conditioned”  set  of  equations 
whose  solution  is  very  unstable  to  small  variations  in  the  known  quantities,  often  leading  to 
unphysical  results.  Successful  inversion  methods  must  therefore  constrain  the  solution  by  various 
assumptions  based  on  physicality.  The  more  expeditious  (though  less  accurate)  technique  utilized  in 
routine  data  processing,  and  the  focus  of  this  report,  employs  the  concept  of  a  geometric  factor 
using  the  approximate  relation  given  by  Eq.  (3)  for  integral  flux.  In  the  following  we  define  a  set  of 
telescope  and  dosimeter  channels  with  geometric  factors  from  which  one  may  compute  approximate 
integral  fluxes. 

A  complicating  factor  in  using  any  of  the  above  algorithms  to  derive  a  flux  spectrum  from  a  set 
of  channel  count  rates  is  that  many  of  the  channels  respond  to  both  electrons  and  protons.  Any 
contribution  to  the  measured  count  rate  by  unwanted  background  will  introduce  an  error  in  the 
inferred  spectrum.  If  one  can  specify  the  background  spectrum  it  is  possible  to  mitigate  this  source 
of  error  by  applying  some  form  of  background  subtraction  algorithm,  although  this  will  not  be 
addressed  in  this  report.  Fortunately,  except  during  the  relatively  rare  time  when  there  is  a  solar 
proton  event  in  progress,  the  outer  zone  belt  is  essentially  free  of  energetic  protons  that  could 
produce  a  background  signal  in  the  “electron  channels.”  Although  the  inner  zone  is  comprised  of 
both  protons  and  electrons,  there  are  some  “proton  channels”  that  have  very  small,  if  any,  electron 
response. 

TR#1  describes  in  detail  the  method  by  which  individual  response  functions  Ri(E)  were 
determined  for  telescope  and  dosimeter  channels  listed  in  Appendix  A.  In  that  report.  Figures  13 
(for  protons)  and  1 5  (for  electrons)  illustrate  that  the  individual  telescope  Ri(E)  are  poor 
approximations  of  either  narrow  Gaussians  or  step  functions,  and  therefore  are  not  appropriate  for 
estimating  either  differential  flux  (through  Eq.  (3))  or  integral  flux  (through  Eq.  (4)).  However,  in 
many  cases  summing  the  response  for  several  channels  yields  a  resultant  R(E)  that  better 
approximates  a  step  function  over  a  limited  energy  range  where  the  expected  particle  flux  is  most 
intense  and  contributes  most  significantly  to  the  integral  in  Eq.  (1).  Such  combined  channels  where 
a  number  of  individual  logic  block  channels  have  been  summed  to  give  a  new  set  of  “standard” 
telescope  channels  T1-T4  are  given  in  Appendix  B.  Appendix  C.l  lists  the  full  set  of  standardized 
flux  channels  (labeled  in  the  first  column)  from  the  dosimeters  (D1-D8)  and  the  telescope  (TOl- 
T09)  forTSXS  (CEASE  sn002).  The  second  column  describes  which  instrument  channel(s)  are 
associated  with  the  standard  channel  label.  The  electron  channel  response  is  defined  by  an  energy 
threshold  Ej  (column  3)  and  geometric  factor  G  (column  4);  the  proton  channel  response  is  defined 
by  Et  (column  5)  and  G  (column  6).  An  analogous  table  exists  for  DSP21  (CEASE  sn007)  in 
Appendix  C.2),  and  DSX  (CEASE  sn004)  in  Appendix  C.3. 
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The  tabulated  geometric  factors  were  determined  using  two  different  methods  referred  to  here  as 
the  Bow  Tie  Analysis  [Van  Allen  et  al,  1974]  and  a  modified  Bow  Tie  Analysis.  Both  methods  use 
the  following  relation  (from  Eqs.  (1)  and  (2),  where  R|(E)  is  the  response  function  for  the  i*  channel 
from  the  set  DO  1 -DOS  orT0l-T09  : 


[R,{E)jXE)dE 

r  mdE 


(5) 


Given  Ri(E),  one  can  solve  for  G,  by  assuming  a  form  of  j(E)  and  an  energy  threshold  Ej.  If  Ri(E) 
were  an  exact  step  function  (with  step  at  Et)  then  this  relation  would  be  an  equivalence  independent 
of  the  spectral  form  of  j(E).  Since  this  is  not  the  case,  its  dependence  on  both  the  spectral  shape  of 
j(E)  and  Ej  introduces  significant  uncertainty  in  the  definition  of  Gj.  To  the  degree  that  the  ramp-up 
of  the  real  response  deviates  from  a  step  function,  the  value  of  G  determined  from  Eq.  (5)  will  be 
sensitive  to  the  chosen  value  of  Ej.  Also,  the  fact  that  the  value  of  Gj  is  dependent  upon  the  spectral 
shape  of  the  measured  flux  distribution  is  obviously  problematic  given  that  the  spectral  shape  may 
vary  over  space  and  time. 

The  Modified  Bow  Tie  Analysis  (MBTA)  is  described  first  since  this  was  the  approach  first 
taken.  The  MBTA  assumes  a  predetermined  range  of  spectral  shape  index  (y)  as  well  as  its  most 
probable  value  and  determines  the  corresponding  values  of  Et  using  a  technique  described  below.  It 
then  computes  G  using  Eq.  (5)  for  each  of  the  three  pairs  of  (y,  Et).  An  uncertainty  in  Gj  may  be 
inferred  from  the  lower  and  upper  limit  of  G  (Gl  and  Gy,  respectively). 

The  range  of  spectral  shape  index  (y)  was  inferred  by  the  differential  flux  spectra  defined  by 
CRRES  instrumentation  (PROTEL  for  protons;  MEA  and  HEEF  for  electrons)  which  had  superior 
energy  and  angular  resolution  relative  to  the  CEASE  instrumentation.  These  results  are  provided  in 
detail  in  Appendix  D.l  (protons)  and  D.2  (electrons),  and  are  summarized  here.  For  both  electrons 
and  protons,  the  CRRES  fluxes  were  fit  to  both  power  law  (j  ~  E'")  and  exponential  (j  ~  exp(-E/Eo) ) 
spectral  shapes,  with  spectral  indices  y  =  n  and  y  =  Eq,  respectively  .  It  was  found  that  throughout 
the  inner  zone  (L  =  1.4-  2),  the  protons  were  best  characterized  by  a  power  law  spectrum  with  n  = 

I  to  4,  with  n  ~  2  at  the  heart  of  the  belt  (L  ~  1 .5).  These  CRRES  results  were  in  relatively  good 
agreement  with  APS  results.  Solar  proton  spectra  may  typically  be  characterized  by  power  laws 
with  n  -  2-3  [Kahler.  2006].  Throughout  the  outer  zone,  electron  spectra  were  found  to  be  best 
characterized  by  an  exponential  spectrum,  although  as  L  ~  6.6  (geosynchronous)  is  approached 
from  lower  L,  power  laws  gave  a  marginally  better  fit.  The  average  characteristic  energy  Eq  is  ~ 

0.35  MeV  for  L  =  3-4  and  decreases  to  ~  0.2  MeV  for  L  =  5-6.  The  average  Eq  for  L  =  3-6  is  ~  0.3 
MeV  with  a  range  of  ~  0.1  to  0.5  MeV. 

The  energy  threshold  Et  for  any  given  R(E)  and  y  is  based  on  the  following  parameter  computed 
as  a  function  of  E  (upper  endpoint  of  the  numerator’s  integral): 

f  Ri,E')j{E')dE' 

S{E,)  =  ^ -  (6) 

I  R{E')j{E’)dE' 


This  is  the  fraction  of  the  total  count  rate  that  is  attributed  to  the  flux  below  a  specified  energy  E. 
The  energy  threshold  Et  is  chosen  to  be  that  energy  for  which  5(E)  first  reaches  (or  exceeds)  the 
threshold  value  ofO.  1.  This  critical  value  of  5  is  chosen  somewhat  arbitrarily,  though  it  provides  a 
quantifiable  criterion  for  defining  Et  that  can  be  consistently  applied  to  all  channel  responses. 
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The  geometric  factor  G,  for  the  i*  channel  is  thus  determined  using  a  specified  spectral  index  to 
define  the  differential  flux  j(E)  in  (5),  and  a  threshold  energy  Et  determined  by  analyzing  the 
function  5(E)  described  above.  Results  using  this  procedure  are  illustrated  in  Figure  1  for  sample 
channels. 


Fig.  1  Results  of  geometric  factor  computation  for  proton  response  to  channel  T5  (left  panels),  and  electron  response  to 
channel  D1  (middle  and  right  panels).  Top  panels  plot  5(E)  vs  E  and  bottom  panels  plot  G  and  R(E)  vs  E.  The  colored 
solid  curves  for  5  (top)  and  the  colored  dashed  curves  for  G  (bottom)  correspond  to  different  spectral  indices  (key  in  top 
panels)  and  threshold  energies  Et  (key  in  bottom  panels).  The  red  curves  correspond  to  relatively  hard  spectral  index 
(low  n,  high  Eo);  the  blue  curves,  to  relatively  soft  spectral  index  (high  n,  low  Eo);  and  the  black  curves,  to  the  more 
typical  spectral  index  as  determined  by  a  survey  of  CRRES  and  NASA  spectra. 


The  channel  T05  proton  response  illustrates  the  sensitivity  of  the  results  to  the  spectral  hardness. 
For  the  power  law  n  =  4  (soft,  steeply  falling  spectrum)  a  substantial  fraction  of  the  counts  come 
from  the  lowest  energy  portion  of  the  response  curve.  Because  the  response  for  T5  is  “two-tiered,” 
the  computed  Et  and  G  for  n  =  1.1  and  2  spectra,  which  are  relatively  flat  and  weight  more  heavily 
the  higher  energy  portion  of  the  response  curve,  are  significantly  larger  than  that  for  the  n  =  4 
spectrum.  The  middle  panels  illustrate  the  electron  response  to  the  D1  channel.  The  response  curve 
(bottom  panel)  shows  a  fairly  ragged  response  below  1  MeV  which  is  attributed  to  a  weak 
bremsstrahlung  signal  with  relatively  poor  statistics  and  does  not  correspond  to  the  penetration  of 
primary  electrons.  For  the  exponential  Eo  =  0.1  MeV  (soft,  steeply  falling  spectrum)  about  80%  of 
the  counts  arise  from  this  part  of  the  response  curve.  This  result  is  compared  to  the  rightmost  panels 
which  show  the  analogous  results  for  D1  where  the  bremsstrahlung  response  is  neglected.  For  very 
soft  spectra  the  uncertainty  in  Et  and  G  can  be  greatly  compounded  by  the  noise  level  near  the 
“tum-on”  point  of  the  response  curves. 

The  deviation  of  some  of  the  response  curves  from  an  ideal  step  function  is  so  great  that  they  are 
worthless  in  terms  of  integral  channels.  Figure  2  illustrates  such  channels.  Channel  D02  has  a 
proton  response  with  such  a  gradual  ramp  that  the  uncertainty  of  Et  is  too  great  to  be  useful. 

Channel  T02  has  a  proton  response  at  ~  1  MeV,  a  broad  minimum  in  response  between  1 .5  and  1 5 
MeV,  and  then  a  ramping  up  to  its  maximum  response.  Channel  T06  has  no  detectable  electron 
response  below  ~5  MeV  and  does  not  flatten  until  beyond  1 0  MeV  which  is  past  the  energy  range  of 
practical  interest.  For  channels  such  as  these,  a  lower  energy  is  listed  in  the  Appendix  C  tables  but 
the  value  of  G  is  labeled  N/A  (Not  Applicable)  because  no  reasonable  G  could  be  computed. 
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Channel  T07  (not  shown)  responds  to  protons  with  very  wide  incident  angles,  missing  the  front 
detector  but  penetrating  the  back  detector.  It  has  become  evident  that  modeling  channel  T07  is 
problematic  since  many  of  the  contributing  protons  penetrate  the  spacecraft  walls  far  from  the 
telescope  cover.  Therefore,  the  range  of  incident  angle  has  been  restricted  to  0°  to  45°  and  channel 
T07  is  ignored.  Ignoring  incident  angles  greater  than  45°  introduces  a  small  error  which  is  relatively 
negligible  compared  to  the  other  uncertainties  discussed. 


/ ;  i 


Figure  2.  Shown  here  (left  to  right)  in  a  format  identical  to  Figure  1  is  the  proton  response  for  channels  D02 
and  T02,  and  the  electron  response  for  T06. 

As  illustrated  by  Figure  3,  the  summation  of  the  dosimeter  LoLET  and  HiLETA  channel  proton 
responses  result  in  a  much  better  approximation  to  a  step  function  then  either  one  by  themselves. 


Figure  3.  Shown  here  (left  to  right)  in  a  format  identical  to  Figure  1  is  the  proton  response  for  channels  D02, 

D04,  and  the  sum  of  D02  plus  D04  (=  DOS). 

A  second  method  that  may  be  used  to  compute  the  geometric  factors  is  the  Bow  Tie  Analysis 
(BTA).  This  method  also  utilizes  equation  (5)  but  is  based  on  a  graphical  technique  for  its  solution. 
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Using  a  predetermined  range  of  the  spectral  index  (y)  and  threshold  energies  (Ej),  several  curves 
parameterized  by  y,  G(Et;  y),  are  computed  and  plotted  as  illustrated  in  Figure  4  where  BTA  is  used 
to  determine  (Ej,  G)  for  electron  dosimeter  channel  D02  assuming  exponential  spectra. 


Figure  4.  Bow  Tie  Analysis  for  TSX  dosimeter  electron  channel  D02.  The  colored  lines  (top  panel)  are 

parameterized  by  y.  For  the  electron  exponential  fits  shown  here,  y  =  Eq  =  0.2  (blue),  0.3  (red),  0.4  (black), 
and  0.5  (green).  The  point  marked  by  the  red  asterisk  in  the  top  panel  corresponds  to  the  threshold  energy  (Et) 
and  geometric  factor  (G)  designated  on  the  response  flinction  plot  (lower  panel). 

The  values  of  Ej  and  G  correspond  (ideally)  to  the  '‘point”  in  the  upper  plot  where  the 
parameterized  lines  converge.  In  practice,  as  Figure  4  illustrates,  the  “point”  of  convergence  is 
generally  quite  spread  out  and  underscores  the  obvious  source  of  uncertainty  in  this  method.  The 
stated  values  have  been  determined  by  an  algorithm  that  determines  the  coordinates  of  all  the  points 
of  intersection  between  the  various  pairs  of  lines.  The  E  coordinates  and  G  coordinates  are  then 
separately  averaged  (and  their  standard  deviation  determined)  to  give  the  solution  (Ej,  G)  which  is 
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then  listed  and  plotted  (with  dashed  lines)  on  the  plot  of  the  response  function  (lower  panel).  One 
has  the  option  of  weighting  the  various  lines  according  to  the  probability  of  encountering  spectra 
with  the  corresponding  spectral  index,  but  this  is  not  done  here.  It  should  also  be  pointed  out  that 
whereas  the  BTA  is  able  to  define  (Et,  G)  for  D02,  the  MBTA  is  not  unless  the  bremsstrahlung 
contribution  at  the  lowest  energies  (<3  MeV)  is  ignored. 

Another  point  to  emphasize  in  comparing  the  results  between  MBTA  and  BTA.  is  that  the  range 
of  solutions  from  BTA  depends  upon  the  range  of  spectral  index  is  used.  In  Figure  5  the  results  for 
the  telescope  proton  channel  T05  are  shown.  The  left  panel,  where  the  4  power  law  indices  used  are 
(1.1  (green),  2.0  (black),  3.0  (red),  and  4.0  (blue)},  shows  a  resulting  Et  =  1 1 .5  MeV.  The  right 
panel,  where  n  =  1.1  has  been  excluded,  shows  a  resulting  Et  =  16.3  MeV.  These  results  may  be 
compared  with  result  from  Figure  1  (left  panel)  where  MBTA  gives  Et  =  22.6.  These  differences 
underscore  the  point  that  although  it  is  more  convenient  to  have  a  specific  energy  and  geometric 
factor  assigned  to  each  individual  channel,  there  are  many  arguments  to  be  made  for  utilizing  the 
full  response  function  and  inferring  a  flux  spectrum  through  some  inversion  algorithm. 


Figure  5.  Comparison  of  Bow  Tie  Analysis  assuming  4  spectral  indices  (left)  and  3  spectral  indices  (right)  for 
telescope  proton  channel  T05. 


The  results  of  both  BTA  and  MBTA  are  summarized  in  Appendix  C.  Whereas  the  results 
illustrated  in  the  Figures  I  to  3  are  based  on  the  earlier  used  upper  limit  energies  for  electrons  (10 
MeV)  and  protons  (200  MeV),  the  final  results  in  Appendix  C  are  based  on  revised  upper  limit 
energies  for  electrons  (50  MeV)  and  protons  (500  MeV)  -  hence,  there  may  be  differences  in  results 
quoted  here  as  compared  to  previous  results. 
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3.  SOFTWARE  DOCUMENTATION 


An  individual  Monte  Carlo  (MC)  run  (MCNPX  or  ACCEPT)  is  made  for  each  pair  of  incident 
particle  energy  (E)  and  angle  (0),  each  run  consisting  of  thousands  of  particle  trajectories  whose 
entry  point  from  the  designated  source  surface  area  is  randomly  chosen.  For  each  of  these  runs,  an 
“input  deck”  is  generated  that  consists  of  the  instrument  geometry  and  material  list  (this  part  is  the 
same  independent  of  the  source  energy/angle  and  is  read  in  as  a  standard  file),  followed  by  several 
lines  which  define  the  source  characteristics  and  the  format  of  the  output  files  (these  lines  are 
generated  by  the  Monte  Carlo  “run”  software).  To  automate  the  process  of  making  hundreds  of 
individual  runs  and  to  facilitate  the  post-MC-run  data  processing  and  bookkeeping  task,  a  number  of 
auxiliary  software  programs  were  developed.  The  various  geometry  files  used  for  creating  the 
“input  decks”  are  summarized  in  Table  1 .  The  specific  programs  for  running  the  Monte  Carlo  codes 
(MCNPX  and  ACCEPT)  and  performing  post-MC  analysis  are  summarized  in  Table  2. 


Table  1.  Input  Deck  Geometry  Files.' 


Sat 

Inst/ 

specie 

Geometry  File  (c:\science\...) 

Description 

Input  to: 

(see  Table  2) 

TSX-5 

Tele  /  e 

ACCEPT\geometry\TEL-TSX  accept  geometry  vOO.txt 

Bare  telescope 

run  accept  tele.for 

ACCEPT\geometry\TEL-TSX  accept  geometry  vOl.txt 

In  frame/case 

DSP21 

Tele  /  e 

ACCEPT\geometry\TEL-DSP  accept  geometry  vOO.txt 

Bare  telescope 

run  accept  tele.for 

ACCEPT\geometry\TEL-DSP  accept  geometry  vOl.txt 

In  frame/case 

DSX 

Tele  /  e 

ACCEPT\geometry\TEL-DSX  accept  geometry  vOO.txt 

Bare  telescope 

run  accept  tele.for 

ACCEPT\geometry\TEL-DSX  accept  geometry  vOl.txt 

Not  used 

All 

Dose  /  e 

ACCEPT\geometry\DD-accept_geometry_vO  1  .txt 

DD1&DD2  in 

case 

run_accept_dose .  for 

TSX-5 

Tele  /  p 

mcnpx\geometry\TEL-TSX  mcnpx  geometry  vOO.txt 

In  case 

run  mcnpx  tele.for 

DSP21 

Tele  /  p 

mcnpx\geometry\TEL-DSP  mcnpx  geometry  vOO.txt 

In  case 

run  mcnpx  tele.for 

DSX 

Tele  /  p 

mcnpx\geometry\TEL-DSX_mcnpx_geometry_vOO.txt 

Tele;DDl;DD2 
in  case 

runmcnpxtele.for 

All 

Dose/  p 

mcnpx\geometry\DOS_mcnpx_geometry_vOO.txt 

DDI;  DD2  in 
case  w/  tele  cyl 

run_mcnpx_dose.for 

1 
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Ta 

j\q  2.  Auxiliary  software  for  Monte  Carlo  analysis. 

Code  Name 

Folder  Location 

(c:\science\...) 

Description 

run_mcnpx_tele.for 

mcnpx\mcnpx_exe\ 

Runs  either  species,  any  satellite,  for 
telescope  -  runs  mcnps.exe,  and  outputs 
p  DSX  E003000  AOO  v03  Ibm.txt,  etc 

runaccepttele.for 

Accept\accept_cease_tele\ 

Runs  electrons,  any  satellite,  for  telescope  - 
runs  accsngl.exe,  outputs 
e  TSX  EOOIOOO  AOO  Ibm.txt,  etc 

responsetele.for 

Response_code\ 

Reads  input  files  (various  versions  and 
formats  depending  who  created  them)  -  for 
either  species  and  any  satellite  and  outputs 
eResponse_DSP_tel.txt, 
eAngularR  DSP  tel.txt,etc 

plot_tele_response.pro 

Response_code\ 

Reads  output  from  response_tele.for  - 
either  species  and  any  satellite,  and  plots 
response  for  channels;  produces  listing  for 
standard  channel  responses 

ru  n_mc  n  pxdose .  for 

mcnpx\mcnpx_exe\ 

Runs  for  either  species,  for  dosimeter-  runs 
mcnpu,  and  outputs 
p_dos_E044000_A00_v0 1  _tbl  .txt,  etc 

Same  output  used  for  any  satellite 

response_mcnpx_dose.for 

Response_code\ 

Reads  in  output  from  run_mcnpx_dose.for, 
creates  LoLET,  HiLETA,  and  HiLETB 
channel  table  for  both  domes 

run_accept_dose.for 

accept\accept  cease_dose\ 

Runs  for  electrons,  for  dosimeter-  runs 
accJ  16.exe,  and  outputs 
e_DOS_E044000_A00_v0  l_out.txt,  etc 

Same  output  used  for  any  satellite 

responseacceptdose.for 

Response_code\ 

Reads  in  output  from  run_accept_dose.for, 
outputs  DD_SN2_eResponse.txt  and 
DD_SN4_eResponse.txt  (all  channels,  both 
dosimeters  -  for  both  TSX5  and  DSP2 1 
(=DSX);  also  separate  files  for  separate 
domes:  DD1(2)  SN2  eResponse.txt,  etc 

plot_dose_response.pro 

Response_code\ 

Reads  dosimeter  output  -  either  species  and 
any  satellite,  and  plots  response  for 
channels 

geometric_factor.pro 

Response_code\ 

Does  modified  Bow  Tie  Analysis;  outputs 
plots  and  tables 

bowtie.pro 

Response  code\ 

Does  Bow  Tie  Analysis;  outputs  plots 
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The  program  run_mcnpx_tele.for  is  designed  to  flexibly  handle  a  variety  of  cases.  The  user 
must  specify  the  CEASE  satellite  (isat),  particle  species  (ispecies),  and  version  number  (iver)  which 
may  vary  depending  upon  (for  example)  the  source  definition  or  executable  version  used.  The  set  of 
energies  and  angles  to  be  run  are  specified  by  an  input  file  sat j9_energies_vxx.txt,  where  sat  is  a  3- 
letter  identifier  for  the  satellite  (TSX,  DSP,  DSXy,  /?  is  a  single-letter  identifier  for  particle  species  {p 
for  proton;  e  for  electron),  and  xx  is  a  2-digit  identifier  for  version  number.  This  file  is  formatted 
with  the  number  of  energies  on  the  first  line,  followed  by  the  set  of  energies;  and  then  the  number  of 
angles,  followed  by  the  set  of  angles.  The  program  loops  through  sets  of  energies  and  angles,  and 
for  each  (E,0)  executes  the  appropriate  code  depending  on  what  flags  are  set.  If  the  imcnpx  flag  is 
set,  the  relevant  instrument  geometry  file  is  read  in  and,  along  with  specific  particle  source 
information,  the  MCNPX  input  deck  is  generated.  MCNPX  is  then  executed  producing  its  standard 
OUTPUT  file.  If  the  iconvert  flag  is  set,  then  the  OUTPUT  file  for  specified  (isat,ispecies,iver,E,0) 
is  read  through  and  only  the  essential  information  (number  of  histories,  incident  energy  and  angle, 
probability  distribution  P(E,0),  average  energy  deposition  per  specified  cell,  etc)  is  extracted  and 
written  to  an  auxiliary  file.  At  this  point  the  P(E,0)  is  specified  for  each  logic  channel  (for  the 
telescope  there  are  9  DFT  x  9  DBT  =  81  logic  channels).  The  imcnpx  and  iconvert  are  independent 
so  that  all  OUTPUT  files  may  be  generated  first  and  then  at  a  later  time  may  be  converted.  The  file 
naming  convention  used  for  the  CEASE  telescope  is  the  following; 

s_SAT_Eeeeeee_Aaa_Vvv_in.txt  -  (INPUT  file) 
s_SAT_Eeeeeee_Aaa_Vvv_out.txt  -  (OUTPUT  file) 
s_SAT_Eeeeeee_Aaa_Vvv_tbl.txt-  (table  of  individual  history  outcomes) 
s_SAT_Eeeeeee_Aaa_Vvv_lbm.txt  -  (summary  table  of  logic  block  matrix  of  P(E,0)  results) 

The  first  character  specifies  the  species  (s  =  p,  proton;  e,  electron);  SAT  species  the  satellite  (SAT  = 
TSX,  DSP,  DSX,  NRL);  eeeeee  specifies  the  energy  E  in  keV;  aa  specifies  the  angle  A  in  degrees; 
and  vv  specifies  the  version  V.  All  numerical  labels  are  padded  with  leading  zeros. 

The  logic  for  run_accept_tele.for  is  similar  to  that  for  run_mcnpx_tele.for,  though  there  are 
some  differences.  In  particular,  the  ACCEPT  codes  require  that  cross-section  tables  be  generated  for 
each  energy,  and  this  is  accomplished  with  run_accept_tele. 

Once  the  logic  channel  P(E,0)  data  are  written,  a  second  program  (response  tele.for)  can  be  run 
which  reads  in  P(E,0)  and  computes  for  the  subset  of  telemetry  channels  the  energy  and  angle 
dependent  effective  area  A(E,0).  Some  telemetry  channels  are  the  individual  logic  block  channels, 
whereas  some  are  summations  of  two  or  more  individual  logic  channels  which  have  been  combined 
because  of  their  similar  response  and  to  conserve  on  telemetry.  Another  file  is  written  that  includes 

T  •  2 

the  energy -dependent  response  R{E)  =  2;r  |  A{E,  6)  sin  ddG  (see  Table  2  for  sample  naming 

0 

conventions). 
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APPENDIX  A  -  PRIMARY  TELESCOPE  AND  DOSIMETER  CHANNELS 


CEASE  /  Telescope 

CEASE  /  Dosimeter 

data  ch# 

channel 

data 

channel 

data 

Channel 

01 

LB(1,0) 

2L 

LB-SUMA 

01 

DDl-LoLET-Flux  (DDILF) 

02 

LB(2,0) 

2? 

LB-SUM  1 

02 

DDl-HiLETA-Flux(DDlHFA) 

03 

LB(3,0) 

23" 

LB-SUM2 

03 

DDl-HiLETB-Flux  (DDIHFB) 

04 

LB(4,0) 

24^ 

LB-SUMB 

04 

DDl-LoLET-DRC  (DDILD) 

05 

LB(5,0) 

25" 

LB-SUM3 

05 

DDl-HiLETA-DRC  (DDIHDA) 

06 

LB(6,0) 

26^ 

LB-SUMC 

06 

DDl-HiLETB-DRC  (DDIHDB) 

07' 

LB(7-8,0) 

Tf 

LB-SUM4 

07 

DD2-LoLET-Flux  (DD2LF) 

08 

LB(0,1) 

LB-SUM5 

08 

DD2-HiLETA-Flux  (DD2HFA) 

09 

LB(1,1) 

29 

LB(5,1) 

09 

DD2-HiLETB-Flux  (DD2HFB) 

10 

LB(2,1) 

30 

LB(5,2) 

10 

DD2-LOLET-DRC  (DD2LD) 

11 

LB(3,1) 

31 

LB(5,3) 

11 

DD2-HiLETA-DRC  (DD2HDA) 

12 

LB(4,1) 

32 

LB(0,2) 

12 

DD2-HiLETB-DRC  (DD2HDB) 

13 

LB(1,2) 

33 

LB(0,3) 

14 

LB(2,2) 

34 

LB(1,4) 

15 

LB(3.2) 

35 

LB(2,4) 

16 

LB(4,2) 

36 

LB(2,5) 

17" 

LB(1,3) 

37" 

LB(6,1) 

18" 

LB(2,3) 

. 38"  ■ 

LB(7,1) 

19 

LB(3,3) 

39" 

LB(8,1) 

20 

LB(4,3) 

40" 

LB(6,2) 

41" 

LB(7,2) 

4P 

LB(8,2) 

iF 

LB(6,3) 

4F 

LB(7,3) 

4F 

LB(8,3) 

Notes 

1.  LB(7-8.0)  =  LB-B32  =  LB(7.0+(8.0) 

2.  CEASE-11  position  17  is  the  sum  of  LB(1,3)  and  LB(2,3);  position  18  is  vacant 

3.  LB-SUMA  =  LB-B46  =  LB  (4.4)+(4,5)+(3,4) 

4.  LB-SUMI  =  LB  (6,1  )+(7.1)+(8.1  H6.2)+(7,2H8,2)+(6.3)+(7,3)+(8,3) 

5.  LB-SUM2  =  LB  (6,4)+(7,4)+(8,4)+(6,5)+(7,5)+(8,5)+(6.6)+(7,6)+(6.7)+(7.7)+(6,8)+(7,8) 

6.  LB-SUMB  =  LB-B49  =  LB  (8.6)+(8,7H8,8) 

7.  LB-SUM3  =  LB  (5.4)+(5.5H5.6H5.7)+(5,8) 

8.  LB-SUMC  =  LB-B51  =  LB (4,6)+(4,7H4.8) 

9.  LB-SUM4  =  LB(3.5)+(3.6)+(3,7>+3,8) 

10.  LB-SUM5  =  LB  (0.4)+(0.5)+(0.6)+(0,7) 

1 1.  Logic  Blocks  numbered  37-45  are  in  CEASE-II  telemetry  only 
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APPENDIX  B  -  TELESCOPE  SUMMED  CHANNEL  DEFINITION 


Primary  Channel 

T1 

T2 

T3 

T4 

LB(1,0)-  1 

X 

LB(2,0)  -  2 

X 

LB(3,0)  -  3 

X 

X 

LB(4,0)-4 

X 

X 

LB(l,l)-9 

X 

X 

LB(2,1)-  10 

X 

X 

LB(3,1)-  1 1 

X 

X 

LB(4,1)-  12 

X 

X 

X 

LB(1,2)-  13 

X 

X 

X 

LB(2,2)-  14 

X 

X 

X 

LB(3,2)-  15 

X 

X 

X 

LB(4,2)-  16 

X 

X 

X 

LB(1,3)-  17 

X 

X 

X 

X 

00 

1 

CN 

X 

X 

X 

X 

LB(3,3)-  19 

X 

X 

X 

X 

LB(4,3)  -  20 

X 

X 

X 

X 
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APPENDIX  C. 


-  CEASE  SN02  (TSX-5)  GEOMETRIC  FACTORS 


Electrons  (Eo  =  0.3;  0.2-0.5  MeV) 

Protons  (r 

1  =  2;  2-4) 

Std 

Ch 

Channel 

Description 

Et  (MeV) 

Ave  (std  dev) 

G  (cm"  sr) 

Ave  (std  dev) 

Et  (MeV) 

G  (cm"  sr) 

DOI 

DDlLF(mBTA) 

1.23 

3.54e-l 

>30 

N/A 

(BTA) 

1.51  (0.07) 

9.18e-l  (1.87e-l) 

>30 

N/A 

D02 

DD2LF 

N/A 

N/A 

>40 

N/A 

2.42  (0.29) 

6.59e-2(5.11e-2) 

>40 

N/A 

D03 

DDIHFA 

1.65 

3.01e-2 

30.7 

1.65e0 

1.88(0.04) 

6.66e-2(8.31e-3) 

26.3  (1.46) 

1 .32e0  (1 .55e-l ) 

D04 

DD2HFA 

2.42 

4.84e-4 

49.3 

1.50e0 

3.16(0.34) 

8.78e-3(7.l6e-3) 

40.4  (3.4) 

1.1  OeO  (2.02e-l) 

D05 

DDIHFB 

6.46 

N/A 

21.6 

6.66e- 1 

6.57(0.01) 

l.45e-3  (2.42e-5) 

16.4(1.9) 

5.l6e-l  (l.29e-l) 

D06 

DD2HFB 

7.41 

N/A 

39.6 

3.30e-l 

7.60  (0.04) 

6.02e-4  (6.02e-5) 

28.9(3.4) 

2.30e-l  (5.86e-2) 

D07 

DD1LF  + 

1.23 

3.65e-l 

32.0 

2.4  leO 

DDIHFA 

1.51  (0.08) 

9.64e-l  (2.00e-l) 

29.5(0.15) 

1.98e0(2.14e-2) 

DOS 

DD2LF  + 

N/A 

N/A ,  5.35e-2 

52.3 

2.60e0 

DD2HFA 

2.43  (0.29) 

6.94e-2  (5.46e-2) 

47.5  (1.27) 

1.99e0(1.12e-l) 

TOl 

16  ch  sum 

0.11 

9.86e-4 

>1 

N/A 

0.14(0.04) 

l.lle-3(1.33e-4) 

T02 

14  ch  sum 

0.15 

7.81e-4 

>1 

N/A 

0.17(0.04) 

8.15e-4(8.62e-5) 

T03 

9  ch  sum 

0.35 

4.70e-4 

81.8 

7.13e-l 

0.42  (0.06) 

6.05e-4  (1. Ole-4) 

71.8(2.44) 

4.66e-l 

T04 

4  ch  sum 

0.57 

2.74e-4 

79.2 

5.74e-l 

0.56  (0.03) 

2.62e-4(2.22e-5) 

67.6  (4.0) 

3.77e-l  (4.72e-2) 

T05 

LBSUM3 

N/A 

N/A 

22.3 

2.42e-2 

16.3(0.6) 

1.68e-2(1.26e-3) 

T06 

LBSUMA 

>4.5 

N/A 

37.5 

1.26e-l 

31.1  (2.4) 

9.68e-2(1.58e-2) 

T08 

LB(3,3) 

0.68 

2.07e-5 

96.6 

4.14e-l 

0.72  (0.06) 

2.38e-5  (3.73e-6) 

77.0(7.2) 

2.41e-l  (4.85e-2) 

T09 

LB(4,3) 

0.83 

6.92e-6 

65.6 

1.46e-l 

0.90(0.10) 

9.l5e-6(2.45e-6) 

52.3  (4.9) 

9.74e-2(1.96e-2) 
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APPENDIX  C.2  -  CEASE  SN07  (DSP21)  GEOMETRIC  FACTORS 


Electrons  (Eo  =  0.3;  0.2-0.5  MeV) 

Protons  (power  n  =  2;  2-4) 

Std 

Ch 

Channel 

Description 

Et  (MeV) 

G  (cm’  sr) 

ET(MeV) 

G  (cm’  sr) 

DOl 

DDlLF(mBTA) 

1.23 

3.61e-l 

>30 

N/A 

(BTA) 

1.51  (0.07) 

9.42e-l  (l.92e-l) 

>30 

N/A 

D02 

DD2LF 

N/A 

N/A 

>40 

N/A 

2.42  (0.29) 

6.76e-2  (5.26e-2) 

>40 

N/A 

D03 

DDIHFA 

1.78 

1.81e-2 

30.3 

I.46e0 

2.02  (0.05) 

3.97e-2  (5.02e-3) 

25.5(1.70) 

I.l6e0(l.66e-I) 

D04 

DD2HFA 

2.82 

5.28e-4 

48.5 

1.27e0 

3.53  (0.26) 

7.77e-3  (4.96e-3) 

38.8  (3.62) 

9.14e-l  (1.84e-l) 

D05 

DDIHFB 

6.46 

N/A 

21.6 

6.66e-l 

6.57(0.01) 

N/A 

16.4(1.89) 

5.l6e-l  (l.29e-l) 

D06 

DD2HFB 

6.54 

N/A 

39.6 

3.30e-l 

6.71  (0.12) 

N/A 

28.9(3.39) 

2.30e-l  (5.86e-2) 

D07 

DD1LF  + 

1.23 

3.65e-l 

32.0 

2.41  eO 

DDIHFA 

1.51  (0.08) 

9.64e-l  (2.00e-l) 

29.5(0.15) 

1.98e0(2.14e-2) 

D08 

DD2LF  + 

N/A 

N/A 

52.3 

2.60e0 

DD2HFA 

2.43  (0.29) 

6.94e-2  (5.46e-2) 

47.5(1.27) 

1.99e0(1.12e-l) 

TOl 

16  ch  sum 

0.1 1 

9.85e-4 

>1 

N/A 

0.14(0.04) 

l.lle-3(1.34e-4) 

T02 

14  ch  sum 

0.15 

7.81e-4 

>1 

N/A 

0.17(0.04) 

8.15e-4(8.65e-5) 

T03 

9  ch  sum 

0.35 

4.16e-4 

65.2 

6.65e-l 

0.37(0.03) 

4.32e-4(3.71e-5) 

58.7  (0.9) 

4.73e-l  (1.49e-2) 

T04 

4  ch  sum 

0.58 

1.59e-4 

60.5 

4.45e-l 

0.56  (0.04) 

1.49e-4(1.38e-5) 

52.1  (3.1) 

3.20e-l  (4.11e-2) 

T05 

LBSUM3 

N/A 

N/A 

22.6 

2.72e-2 

16.7  (0.6) 

1.90e-2(l.50e-3) 

T06 

LBSUMA 

>4.5 

N/A 

35.1 

6.29e-2 

28.5(2.3) 

4.80e-2  (8.39e-4) 

T08 

LB(3,3) 

0.66 

1.21e-5 

92.6 

3.27e-l 

0.70  (0.07) 

1.36e-5(2.43e-6) 

71.9(7.4) 

1.90e-l  (4.23e-2) 

T09 

LB(4,3) 

0.83 

4.40e-6 

53.5 

2.08e-l 

0.93(0.10) 

6.45e-6(1.85e-6) 

43.7  (3.9) 

1.50e-l  (2.85e-2) 
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APPENDIX  C.3  -  CEASE  SN04  (DSX)  GEOMETRIC  FACTORS 


Electrons  (Eo  =  0.3;  0.2-0.5  MeV) 

Protons  (power  n  =  2;  2-4) 

Std 

Ch 

Channel 

Description 

ET(MeV) 

G  (cm*  sr) 

Et  (MeV) 

G  (cm“sr) 

DOl 

DDILF(mBTA) 

1.23 

3.61e-l 

>30 

N/A 

(BTA) 

1.51  (0.07) 

9.42e-l  (l.92e-l) 

>30 

N/A 

D02 

DD2LF 

N/A 

N/A 

>40 

N/A 

2.42  (0.29) 

6.76e-2(5.26e-2) 

>40 

N/A 

D03 

DDIHFA 

1.78 

1.81e-2 

30.3 

I.46e0 

2.02  (0.05) 

3.97e-2  (5.02e-3) 

25.5  (1.70) 

I.l6e0(l.66e-1) 

D04 

DD2HFA 

2.82 

5.28e-4 

48.5 

I.27e0 

3.53  (0.26) 

7.77e-3  (4.96e-3) 

38.8(3.62) 

9.14e-l  (1.84e-l) 

DOS 

DDIHFB 

6.46 

N/A 

21.6 

6.66e- 1 

6.57(0.01) 

N/A 

16.4(1.89) 

5.16e-l  (1.29e-l) 

D06 

DD2HFB 

6.54 

N/A 

39.6 

3.30e-l 

6.71  (0.12) 

N/A 

28.9  (3.39) 

2.30e-l  (5.86e-2) 

D07 

DD1LF  + 

1.23 

3.65e-l 

32.0 

2.4  leO 

DDI HFA 

1.51  (0.08) 

9.64e-l  (2.00e-l) 

29.5  (0.15) 

1.98e0(2.14e-2) 

DOS 

DD2LF  + 

N/A 

N/A 

52.3 

2.60e0 

DD2HFA 

2.43  (0.29) 

6.94e-2  (5.46e-2) 

47.5  (1.27) 

1.99e0(1.12e-l) 

TOI 

16  ch  sum 

0.11 

1.22e-3 

>1 

N/A 

0.13  (0.01) 

1.29e-3(4.69e-5) 

T02 

14  ch  sum 

0.16 

9.57e-4 

>1 

N/A 

0.17(0.01) 

9.77e-4(3.60e-5) 

T03 

9  ch  sum 

0.37 

5.44e-4 

78.6 

6.l5e-l 

0.35(0.01) 

5.07e-4(7.l7e-6) 

69.1  (1.7) 

4.06e-l  (2.l6e-2) 

T04 

4  ch  sum 

0.59 

l.94e-4 

72.2 

3.89e-l 

0.53  (0.01) 

l.58e-4(6.36e-6) 

60.7(4.1) 

2.61e-l  (3.75e-2) 

T05 

LBSUM3 

N/A 

N/A 

25.3 

l.20e-2 

19.0(1.2) 

7.90e-3(l.03e-3) 

T06 

LBSUMA 

>9 

N/A 

42.6 

4.38e-2 

34.4  (2.4) 

3.20e-2  (4.70e-3) 

T08 

LB(3,3) 

0.68 

1 .28e-5 

101.8 

2.84e-l 

0.64(0.01) 

1.09e-5  (4.15e*7) 

78.7(8.2) 

1.58e-l  (3.56e-2) 

T09 

LB(4,3) 

0.87 

4.80e-6 

62.9 

1.64e-l 

0.84(0.01) 

4.27e-6(1.48e-7) 

50.8  (4.6) 

I.l2e-I  (2.l8e-2) 
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APPENDIX  D.l  -  CRRES  SPECTRAL  ANALYSIS  (PROTONS) 
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Fig.  6  -  Shown  are  proton  spectra  from  APSmax  (left)  and  CRRESPRO-Quiet  (right)  for  L  -  1 .5.  The  spectra  are  fit  to  a  power  law 
over  the  range  10  to  100  MeV  and  the  fitted  power  law  index  (n)  is  provided  in  the  legend. 
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Fig.  7  -  Shown  are  proton  spectra  from  APSmax  (left)  and  CRRESPRO-Quiet  (right)  for  L  -  2.  The  spectra  are  fit  to  a  power  law  over 
the  range  10  to  100  MeV  and  the  fitted  power  law  index  (n)  is  provided  in  the  legend. 
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APSmax  AP8max 


Fig.  8  -  Shown  are  the  computed  proton  (from  APSmax  model)  power  law  index  (left)  and  exponential  index  (right)  plotted  versus  L. 


Fig.  9  -  Shown  are  the  computed  proton  (from  CRRESpro-Quiet  model)  power  law  index  (left)  and  exponential  index  (right)  plotted 
versus  L. 
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APPENDIX  D.2  -  CRRES  SPECTRAL  ANALYSIS  (ELECTRONS) 


Fig.  10  -  Shown  are  CRRES  electron  spectra  from  orbit  72,  L  -  4.0  (left)  and  orbit  88,  L  -  3.0  (right).  The  spectra  are  fit  to  a  power 
law  (dotted  line)  and  exponential  (dashed  line),  with  the  fitted  power  law  (n)  and  exponential  (Eo)  indices  provided  in  the  legend. 
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2  5  +  5  6  7  8 


Power  Law  Index  (n) 


Fig  11.  -  Histograms  showing  the  frequency  of  occurrence  for  a  range  of  power  law  indices  (top)  and 
exponential  indices  (bottom)  computed  from  CRRES  electron  spectra,  for  three  different  L  ranges  (shown  in 
legend). 
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